Striped Multiferroic Phase in Double-Exchange Model for Quarter-Doped Manganites 
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The phase diagram of quarter-hole-doped perovskite manganites is investigated using the double- 
exchange model. An exotic striped type-II multiferroic phase, where 25% of the nearest-neighbor 
spin couplings are orthogonal to each other, is found in the narrow-bandwidth region. Comparing 
with the spiral-spin ordering phase of undoped manganites, the multiferroic Curie temperature of 
the new phase is estimated to be ~ 4 times higher, while the ferroelectric polarization is similar 
in magnitude. Our study provides a path for non-coUinear spin multiferroics based on electronic 
self-organization, different from the traditional approach based on superexchange frustration. 

PACS numbers: 75.80.-|-q, 75.47.Lx, 71.10.Hf, 75.30.Kz 



Introduction. Multiferroics, materials where ferroelec- 
tric (FE) and magnetic orders coexist, have attracted 
enormous interest due to their technological relevance 
and fundamental science challenges 1]. Based on the 
microscopic origin of the FE polarization {P), the mul- 
tiferroics can be classified into two families Type-I 
multiferroics, where ferroelectricity and magnetism have 
different origins, are often good ferroelectrics with high 
FE transition temperatures and large P. However, the 
coupling between magnetism and ferroelectricity is usu- 
ally weak. In contrast, in the type-II multiferroics (mag- 
netic multiferroics) , such as iJMnOs (i?=Tb and Dy) [3| , 
the ferroelectricity is caused by, and thus it is strongly 
coupled with, a particular magnetic order. The pursuit 
for type-II multiferroics with higher critical temperatures 
Tq and larger P is one of the most fundamental chal- 
lenges of modern condensed-matter physics [1] . 

The manganite family provides fertile ground for both 
type-I (e.g. BiMnOs and hexagonal YMnOa) and 
type-II (e.g. .RMnOa perovskites, i?=Tb, Dy, Ho, or 
Eui_a;Ya;) multiferroics. Recent theoretical studies pre- 
dicted additional typc-I multiferroics in manganites, e.g. 
Pro.6Cao.4Mn03, CaMnOa, and BaMnOa [5]. However, 
since among type-II multiferroics the manganite multifer- 
roics show the largest P (e.g. ~ 0.6-2.5 x 10'^ /iC/m^ for 
TbMnOa and DyMnOg [3], - 10'' ^iC/m^ for HoMnOg 
0), it would be even more exciting if new type-II man- 
ganite multiferroics with higher Tq could be discovered. 

Although the manganite phase diagrams have been 
extensively studied [7|, exotic new phases may still be 
waiting to be discovered, particularly in the narrow- 
bandwidth limit where the multiferroic manganites are 
located. Therefore, it is important to continue the de- 
tailed investigation of this narrow-bandwidth regime. 
Along this direction, recent studies addressed the sys- 
tematic phase diagram for undoped manganites using the 
two-orbital double-exchange (DE) model [1| . 

Compared with the undoped case, the doped multifer- 



roics, which have not been widely studied, can provide a 
fruitful playground since one extra degree of freedom, the 
charge, is now active. It is already known that doping- 
driven electronic self-organization can be very important 
to the emergence of hi^- Tc superconductivity and colos- 
sal magnetoresistance [9|. Thus, it is natural to explore a 
similar path in multiferroics. For instance, very recently, 
the Ca-do ped BiFeOs and TbMnOa were studied exper- 
imentally [lol . flU . In this Letter, theoretical investiga- 
tions are extended to the doped manganites, revealing 
a type-II multiferroic phase at quarter-doping that was 
missing in previous theoretical studies jl2l |. 

Model and methods. In this letter, the phase diagram 
of quarter-doped manganites will be studied using the 
two-orbital DE model. The Hamiltonian is: 



H — HuE{to) + Hse{Jaf) + -fflattico(A), 



(1) 



which includes the two-orbital DE hopping term -ffoE, 
the nearest-neighbor (NN) superexchange (SE) energy 
HsEi and the lattice contribution iJiatticc that contains 
the electron-phonon coupling and lattice's elastic energy. 
Details of this well-known Hamiltonian can be found in 
previous publications 0, [1] . The energy unit will be the 
DE hopping to (~ 0.3-0.5 eV). Besides to, there are two 
main parameters in this Hamiltonian: Jaf for the SE, 
and A for the electron-phonon coupling. 

Three numerical methods have been employed to cross 
check the results. The first technique is the zero tem- 
perature (T) variational method where the ground-state 
energies of several phases are compared. This calculation 
can be applied on both finite- and infinite-size lattices. 
However, a set of candidate phases must be preselected. 
To avoid such a possible bias, the finite-T Monte Carlo 
(MC) simulation is also here performed, on finite-size 
clusters, to confirm the phase diagram obtained by the 
variational method. The typical clusters, with periodic 
boundary conditions, used here are 2-dimensional (2D) 
L X L {L — 8 and 16). Some hidden phases or phase- 
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FIG. 1: (Color online) The zero-T (2D lattice) phase dia- 
gram of quarter-doped manganites. SOS denotes the novel 
multiferroic phase. Full dots denote the couplings where vari- 
ational results were confirmed by the MC simulation J^J and 
the zero-T relaxation using a L = 8 cluster. Open dots are 
cases where the MC simulation did not provide a clear answer 
due to strong competition of metastable states. 

separation tendencies missed by the first method can be 
revealed in the MC simulation. Finally, the zero-T re- 
laxation technique is used to optimize the ground state 
spins' and lattice' patterns on a finite-size cluster. The 
second and third methods are only applied to 2D lattices 
in the current effort because they are computationally 
intense. The three processes are cycled several times to 
improve the accuracy of the phase diagram. 

Results and Discussions. The 2D ground state phase 
diagram of the quarter-doped DE model is shown in 
Fig. 1 [14]. A prominent FM metallic orbitally-disordered 
phase occupies the small Jaf region, in agreement with 
the phase diagram of large and intermediate bandwidth 
manganites. However, our most important result is re- 
vealed when Jaf is increased, keeping A small: in this 
regime a new multiferroic phase (dubbed "SOS" as ex- 
plained below) is found. Its existence is confirmed by all 
the three methods used here. In addition, the C1/4E3/4 
phase, predicted by Hotta et at p^H, appears when A~l. 
Finally, a G-AFM region is found when both A and Jaf 
are large enough [l^ . Near the boundaries between these 
robust phases the spin patterns are difficult to identify 
and phase separation may exist there (shaded regions in 
Fig. 1) fl6i|. These regions are beyond our current com- 
putational capabilities. Thus, here the emphasis will be 
only on the clearly identified states in the phase diagram. 

Let us now analyze the physical properties of the SOS 
phase, comparing properties with those of the C1/4E3/4 
phase. The spin patterns of these two phases are shown 
in Fig. 2(a-b) (both of them are perfectly AFM along 
the c-axis). For the SOS phase, the most prominent 
features are the diagonal stripes where spin domains 
with orthogonal orientations merge. For this reason, 
this phase is called "Spin-Orthogonal Stripe" (SOS). 
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FIG. 2: (Color online) Spin patterns of the (a) SOS and (b) 
C1/4E3/4 phases. The magnetic unit cell of the SOS phase 
(red square) is 4\/2 x 4\/2 (4 x 4 in the orthorhombic nota- 
tion). Shown also are the two possible local P from the DM 
mechanism (red arrows). The total P is, thus, perpendicular 
to the stripes [TJ- (c) The (common) spin structure factor 
of these two phases. The characteristic peaks appear at (1/4, 
1/4), (1/8, 3/8), and (3/8, 1/8) and their equivalent positions 
in the q + (1/2, 1/2) region. The ideal amplitudes of these 
three peaks are in the ratio of 2 : 1 : 1. Density-of-states of 
the (d) SOS and (e) C1/4E3/4 phases (at A = 0, fixing the 
spin pattern). Both of these phases have a gap (0.6, ~ 0.2-0.3 
eV) at the Fermi level. 

This spin pattern gives rise to a FE P at the stripe 
boundaries via the Dzyaloshinskii-Moriya (DM) interac- 
tion {Pi, {Si X Sj), where eij is the unit vector 
connecting the NN spins 5^ and Sj ) [l^l . The two types 
of local P that appear in the SOS state are shown in 
Fig. 2(a). The total P arises from the sum of these local 
P, and it is perpendicular to the stripes. The local \P\ 
at the stripe boundaries should be larger than those in 
TbMnOa and DyMnOa, since here Si x Sj is the largest 
possible in noncoUinear spin orders. However, consid- 
ering that only 25% NN spins are orthogonal, then the 
global \P\ of the state should be similar {^ 10"^ ^C/m^) 
as in TbMnOs and DyMnOa. In contrast, the C1/4E3/4 
phase has all spins coUinear, forming FM zig-zag chains 
with alternative C- and E-AFA4 sections [19|. 

Note that the spin structure factors {S{q) = -5^ • 

S'i+pexp[z27r(f • f]) of the SOS and C1/4E3/4 phases are 
exactly the same (Fig. 2(c)) [20]. The reason is that the 
spin pattern of the SOS phase can be viewed as the com- 
bination of two C1/4E3/4 patterns: by decomposing the 
spin vectors in Fig. 2(a) along the x and y axes, both the 
X and y components by themselves form the C1/4E3/4 
state, with just a relative shift between the two C1/4E3/4 
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FIG. 3: (Color online) Finite-T MC results (Jafm=0.18, 
A=0) for the SOS phase, on the L=8 cluster. Left axis: P 
per site is calculated using Yl<i ^ (Si x Sj)]/L^, phe- 

nomenologically. Thus, \P\ per site at zero-T is 72/4. The 
MC average is done on |-P|^, to avoid strong finite-size fluc- 
tuations in P, and (|-P|^)^^^ per site is shown. Right axis: 
Amplitudes of characteristic peaks in the spin structure factor 
and their summation are normalized to their zero-T values. 
Dashed lines indicate the size-dependent high-T asymptotic 
values of spin structure factors and \P\ (i.e. background). 
Real-space spin correlations at the maximum length C{L/2, 
L/2) are also shown. Inset: same properties but on L—16 
lattices. 

patterns formed by this procedure. Therefore, care must 
be taken in separating the SOS and C1/4E3/4 states ex- 
perimentally because they are intimately connected. 

It is also important to confirm that the SOS phase is an 
insulator, as required by ferroelectricity. The density of 
states on the infinite-size lattice indeed shows that both 
the SOS and C1/4E3/4 phases have a gap at the Fermi 
level (Fig. 2(d-e)). Thus, it is safe to predict that the 
SOS state is a type-II multiferroic. 

The theoretical prediction of the SOS phase is excit- 
ing for its new mechanism for noncoUinear spin ordering, 
which can improve the Tq. Traditionally, noncoUinear 
spirals are formed by the competition between the NN 
and next-nearest-neighbor (NNN) SEs, namely employ- 
ing the SE frustration [8j,l21)]. However, because the NNN 
SE is usually weak, it becomes difficult to improve Tq 
based on such a frustration mechanism. Another known 
path for noncoUinear spin order is the DM interaction 
[18|, but it is also very weak and thus can only induce 
small spin angles. In contrast, the 90° spin angle in our 
SOS phase is not driven by the NNN exchange frustration 
nor the DM interaction. Instead, the noncoUinear spin 
structure of the SOS phase is driven by the competition 
between DE and NN SE. While in cuprate and nicke- 
late stripes the spins across the stripes are AFM ordered 
(TT-shifted) , in manganites there are also FM tendencies 
caused by DE and, thus, a reasonable compromise at the 
stripes is to have 90° spin angles. Since both DE and 
SE are NN interactions, much stronger than NNN ones, 
a higher SOS Tq is expected. 
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FIG. 4: (Color online) Zero-T charge/orbital ordering of the 
SOS and C1/4E3/4 phases. The unit cell of the charge/orbital 
ordering block shown here is only 1/4 of that of the spin or- 
dering (dotted regions in Fig. 2). The radius of the circles 
are proportional to the local eg charge density. The high (H) , 
medium (M), and low (L) charge density regions are parti- 
tioned by the broken lines. For the SOS phase (A = 0), 
the Cg density in the H (L) region is 5/6 (2/3). For the 
C1/4E3/4 phase at A = 0.8, the eg density in the H/M/L 
sites is 0.95/0.73/0.59, respectively. 



To confirm this expectation, the finite-T MC results, 
including the EE \P\, the spin structure factors, and the 
real-space spin correlation at the maximum distance, are 
shown in Fig. 3. All these results give a very similar 
Tq , here defined as the T where the magnetic correlation 
length is as large as the clusters studied. For the L=8 lat- 
tice, the Tc of the SOS phase is about 0.025-0.030, which 
is ~3.85-4.6 times the Tc of the spiral state of previous 
investigations (0.0065, obtained using the same MC tech- 
nique on a L=12 lattice [s^]). Even for the L=16 lattice, 
the SOS Tc is about 0.0225-0.025, stiU much higher than 
in the spiral state and showing that size effects are small. 
Therefore, the real Tc of the SOS phase is estimated to 
be ~ 100 K (4 times the Tc of typical spiral mangan- 
ites). It is reasonable to expect the discovery of other 
multiferroic noncoUinear phases based on this electronic 
self-organization mechanism with even higher Tc's. 

Charge/Orbital order and Coulomb interactions. Both 
the SOS and the C1/4E3/4 phases have a complex 
charge/orbital order, forming stripes along the 6-axis (see 
Fig. 4). For the SOS phase, holes accumulate at the 
diagonal stripe boundaries where P originates. For the 
C1/4E3/4 phase, electrons accumulate at the E-AFM cor- 
ner sites while holes accumulate at the C-AFM bridge 
sites 12|, in agreement with the pure E-AFM (undoped 
case) and C-AFM (over half-doped cases) properties. 
Also, note that the charge disproportion in C1/4E3/4 is 
more prominent than that in the SOS phase. Thus, if the 
Coulomb interaction is considered, the C1/4E3/4 phase 
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will be further suppressed in the phase diagram. Using 
the mean- field approximation [23 |. it can be shown that 
the phase boundary SOS-C1/4E3/4 can be shifted from 
A « 0.7 to 1 by incorporating the on-site Hubbard inter- 
action. The orthorhombic distortion (lattice shrinkage 
along the c-axis and expansion in the a-b plane) will also 
shift the phase boundary SOS-C1/4E3/4 toward A '--^ 1. 

The challenges for the experimental realization of the 
SOS phase are substantial and will require much care: 
(1) The perovskite manganites must be synthesized in 
the narrow-bandwidth limit. To pursue such a large 
JAFu/to (> 0.17), crystal with small lattice constants 
are needed (smaller than in Pr3/4Cai/4Mn03 [13] )■ (2) 
Another challenge is to avoid quenched disorder that may 
transform long-range ordered phases into glassy states 



111 |2j|. These challenges may be solved by carefully 



choosing divalent doping cations (not restricted to alka- 
line earths) similar in size to those of small rare earths 
(e.g. Tb, Dy, and Ho). Also, recent developments on 
complex oxide heterostrucutures may be helpful by cre- 
ating A-site ordered manganites [2HI and controlling the 
crystal lattice using strain/stress. 

Conclusions. Here, a striped type-II multiferroic phase 
of the two-orbital double-exchange model for quarter- 
doped manganites has been reported. This phase is 
spin/charge/orbital ordered and has ferroelectricity due 
to the presence of 90° NN spin angles at spontaneously 
formed diagonal stripes. This noncoUinear spin-ordered 
state potentially induces a multiferroic Tc higher than 
those of spiral-spin ordered undoped manganites. 
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